Incoherent Mollow triplet 
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A counterpart of the Mollow triplet (luminescence lineshape of a two-level system under coherent 
excitation) is obtained for the case of incoherent excitation in a cavity. Its analytical expression, in 
excellent agreement with numerical results, pinpoints analogies and differences between the conven- 
tional resonance fluorescence spectrum and its cavity QED analogue under incoherent excitation. 
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Mollow [T] discovered a striking type of spectral shape 
in the resonance fluorescence problem, where an atom 
is irradiated by a strong laser beam. The celebrated 
Mollow triplet [2], that results from transitions between 
atomic states that are dressed by the coherent light 
field, has since been a testbed of nonlinear optics. It 
stands as one of the fundamental spectral shapes of light- 
matter interaction, maybe second only to the Rabi dou- 
blet. Although the Mollow triplet is rooted in quan- 
tum physics and bears much quantum features itself, 
it arises from a fully classical light field. Its Hamilto- 
nian, in the rotating frame of the laser and at resonance, 
simply reads Hr, — f2z,(er + cr^), with £l\ the laser in- 
tensity and a the only quantum operator, namely, the 
two-level system annihilation operator. Including the 
spontaneous decay of the emitter, in the Lindblad form 
£<r(p) = (2a pa^ — a^ap — per' a), leads to a master equa- 
tion dtp = i[p, Hr] + ^C a (p) from which one obtains the 
famous Mollow triplet lineshape: 
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It is composed of an elastic scattering peak, the Dirac 
S function, and the triplet itself, with a central Lorentzian 
peak of full width at half maximum (FWHM) -f a and two 
satellite peaks at the symmetric positions ±^R(Rt,) with 
Mollow splitting R L = v / (2fi L ) 2 - (7 ff /4) 2 and FWHMs 
37 CT /2. This structure was observed a long time ago with 
atoms [3| and more recently also in a variety of solid 
state systems [IHE], with, as befits the above description, 
coherent excitation. 

In this text we consider a close counterpart of this fun- 
damental system, where the light field is initially fully 
quantized, and becomes continuous as a result of an in- 
coherent and continuous pumping that feeds the system 
with a very large number of photons. This situation is 
realized — as for quantization of the light field — in cav- 
ity QED [5], where quanta of a trapped standing waves 
(the photons) can be brought to interact with an isolated 
emitter, and — as for the incoherent pumping — with semi- 
conductor microcavities |10j . where excitations are con- 



tinuously poured into the system with no external coher- 
ence fed in by a driving field. The role of the emitter is, in 
this case, played by a quantum dot placed in the antinode 
of the microcavity field. In the cavity QED version of the 
resonance fluorescence physics, the system is described by 
the Jaynes-Cummings Hamiltonian (still at resonance): 
H = g(a) a + aa^) with the cavity mode also quantized 
through the boson operator a. Cavity decay j a as well 
as incoherent pumping P a are described like before with 
a master equation d t p — i[p, H] + ' 1 fC a (p) + ^f-C a t(p) 
where p is now the density matrix for the combined two- 
level-emitter/cavity system. We assume that -f a <c P a , 
so much so that we can neglect it for simplicity (it is al- 
ready very small in a typical system and we consider the 
case of large pumping). We also assume the condition 
of strong-coupling, where g 3> 7 Q /4. This regime was 
recently observed in the cavity emission as a transition 
from the Rabi doublet to a single lasing line [IT]. Since 
this transition has been claimed remaining in strong- 
coupling, it results from climbing the Jaynes-Cummings 
ladder il2:, and as such, is a successful realization of 
quantum nonlinearities in these systems. The impor- 
tance of this breakthrough for microcavity QED is how- 
ever hindered by such a simple manifestation, in partic- 
ular since other mechanisms can also result in a simi- 
lar behaviour of Rabi splitting collapse without entering 
the quantum nonlinear regime |13) . Quantum features 
are generally better observed when probing the quantum 
emitter, rather than the cavity, whose close connections 
with the classical oscillator tend to surface rapidly and 
dominate strongly. The theoretical description, which is 
straightforward in the low excitation regime even when 
solving the system exactly [T3], becomes computation- 
ally demanding when the lasing regime is approached, 
but good approximations can be sought [2]. In this 
text, we consider this physics of the highly nonlinear 
regime of the microcavity-QED system, that is, Jaynes- 
Cummings physics under a strong incoherent pumping. 
We find that, in good systems by the standard of today, 
a new type of Mollow triplet is obtained in the direct 
QD emission spectrum. It is a close counterpart of the 
classical Mollow triplet where light is described by a clas- 
sical field [T], whereas it is here described as numerous 
quanta of the cavity mode. The coherence is acquired 
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FIG. 1: (Color online) Exact numerical solution (points) 
and their analytical approximation, Eqs. u3p , (lines), 
for n a (blue/circles), n„ (brown/squares) and g^ 2 ' (pur- 
ple/triangles), for 7 a = O.lgr, as a function of pumping P a /g. 
The analytical solutions become unphysical when P a — k ct 
(here at 40<?), where n a = 0, n„ — 1 and g^ 2 ' diverges. They 
are very good approximations in the region of interest, Eq. (Mk. 



through the strong-coupling with the dot, resulting in 
striking variations from the case where it is provided by 
an external laser. We now describe them analytically. 

Mollow regime - - Whereas only one parameter (in- 
tensity) fully describes the light in Mollow's descrip- 
tion, the Jaynes-Cummings picture requires from the 
start to take into account an infinite number of cor- 
relators between the fields, that we can however re- 
late to each other 12 : (a^ a n ~ l a) = ij^loy a n ) and 



low pump, but still enough to be beyond the quantum 
regime [12], i.e., j a < Pa <C K a , the cavity popula- 
tion increases linearly with pumping, with a half occu- 
pied dot. This is the most effective region for accumu- 
lation of photons in the cavity (the so-called one-atom 
laser 15 ), with little disruption from incoherent pro- 
cesses. At smaller efficiency F a , the dot occupation also 
increases linearly with pumping, quenching the linear in- 
crease of the cavity population. F a represents, there- 
fore, the degree to which the dot pumping succeeds in 
populating the dot itself, against the coherent exchange 
of population that feeds the cavity with efficiency F a . 
These expressions are thus valid until the dot population 
is fully inverted, at P max ~ K&, then, the self-quenching 
dominates the dynamics, emptying the cavity that goes 
to a thermal state. The maximum population of the cav- 
ity, max(n a ) ~ g 2 /(2j 2 ), is reached at the intermediate 
rate P a sw K a /2. This identifies the regime of interest for 
the observation of the Mollow triplet, where the cavity 
field is intense (n a S> 1) and coherent (with a Poissonian 
photon distribution, T[n] = e~" a n"/n! and <?' 2 ' = 1): 



7 Q < g < Pa- < K a 



(4) 
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Now that we have a good and analytical description of 
the populations, we turn to the optical emission spec- 
trum, that we show can be obtained in equally good 
approximations. The dot emission reads n a irSi nc {uj) 



(TV) = [P CT ( a t n a n *) - 7 a VV)]/[P CT + 3? io°° (ct 1 (0)cr(T))e iLU ^d-r. We compute the two-time cor- 



7 a (n — 1)] (all others being zero in the steady state) 
From this follows a first relation for the populations 
of the modes, n a — {era) and n a — (a* a), namely 
n a = (P a — 7 a n a )/ P a . This also allows to obtain a self- 
contained equation for (at a n ): 
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where K a = Ag 2 /j a is the Purcell rate of transfer of pop- 
ulation from the dot to the cavity mode. This recurrence 
equation allows to compute (<v a n ) for all n as a func- 
tion of n a only. The solution for n = gives a good 
approximation for the region where the cavity field be- 
haves classically: 
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The quality of this approximation is seen in Fig. [T] where 
it is compared with the exact solution, computed nu- 
merically [T5]. The second order coherence function g^ 
also admits a closed-form expression (not given here but 
plotted in Fig. fl]) which is unity in good approximation. 
The expressions Eqs. (pi for the populations have a clear 
physical meaning, that can be discussed in terms of two 
parameters: the "cavity feeding", F a — P a -/(2j a ), and 
the " dot feeding" , F a = (P a — 7a)/ft<ri efficiencies. At 



relator (crt(o)cr(r)) in two steps: first, we solve the master 
equation in the steady state, finding the density matrix 
elements p m ,i;n,j (for m, n e N and i, j £ {0, 1}, pho- 
ton and exciton indexes, respectively). For the range of 
parameters of interest, we show that they can be analyt- 
ically expressed only in terms of the photon distribution 
T[n\. Second, we apply the quantum regression formula. 
1. Steady state density matrix -- We consider only 
elements that are nonzero in the steady state: the popu- 
lations Pi[n] = p n ,i-n,i with i = 0, 1, corresponding to the 
probability to have n photons with (j>i) or without (po) 
exciton, and the off-diagonal terms qi[n] — ^(p n ,0;n-i,t), 
corresponding to the coherence between the states \n, 0) 
and \n — 1,1). The master equation now reads: 

d t po[n + 1] = T> phot {pr)[n + 1]} (5a) 

- P a p [n + 1] - 2gy/n + l qi [n + 1] , 
d t pi[n] = T> P hot{pi[n]} (5b) 

+ P a {T[n] -pi[n\) + 2gVn + l qi [n + 1] , 

dtQi[n + 1] = ~Dphot{<li[n + 1 ]} (5c) 

p 

- ^f®t n + !] + 9vn + l(po[n + 1 }- PiN) , 

where we have separated the photonic dynamics into a 
superoperator 2? p hot- Given that it is much slower than 
the dot dynamics, one can solve the steady state ignoring 
£>phot [U]. The photon distribution, T[n] = Po[n]+Pi[n], 



remains unperturbed during the excitation and inter- 
action with the dot. Equations (pi) then admit solu- 
tions in terms of Tin], i.e., Po[n] ~ p "y" , ', u T[n] and 

*N ~ P ^ 9 2 ^"„ T[n-l] where k q = 4g 2 /-P<r is the Purcell 
rate of transfer of population from the cavity to the dot. 
Our approximations of large intensities imply n + 1 « n. 
,2. Two-time correlator and spectra — The two-time 
correlator can be expressed as a sum (a^ (O)cr(r)) = 
Y^=o Q[ n ]( T ) i wnere Q[ n ] an d other functions So,i[n] 
and V[n] are defined through the quantum regression for- 
mula by coupled differential equations (n > 0): 



d T Q[n] = T> p hot{Q[n]} 
P 



(6a) 



[n] + ig(y/nSi[n] - \Jn + lS [n + 1]) , 



d T S [n + 1] = V phot {S [n + 1]} (6b) 

- P a S Q [n + 1] + ig{y/nV[n + 1] - y/n + lQ[n]) , 
drSi[n]=V phot {S l [n}} (6c) 

+ P a (X[n] - Si[n]) - ig(Vn + lV[n + 1] - ^JnQ[n]) , 
d T V[n + l]=V phot {V[n+l]} (6d) 

- -fV[n + 1] + *s(%/nS [n + 1] - Vn + lSx[n}) . 

They are, like for the single-time dynamics, separated 
into a slow photonic dynamics embedded in a super- 
operator Dphot that is r-independent in good approxi- 
mation, and a fast exciton and coupling dynamics that 
we can solve analytically. Moreover, we have introduced 
the steady state function X[n] = S [n](0) + Si[n](0), in 
analogy with T[n]. The initial conditions in Eq. (J6J), 
are the steady state values So[n + 1](0) = iqi[n + 1], 
Si[n](0) = 0, Q[n](0) = p 1 [n] and V[n + 1}{0) = (there- 
fore, X[n] = ig;[n]). After some long, but straightforward 
algebra, we can find the expression for Q[?i](t) in terms 
of Po.iN an d q_i[n], which, in turn, are expressed in terms 
of the statistics T[n] (as shown previously). This allows 
to compute a closed-form solution for (ct^(O)ct(t)), which 
is however lengthy and not worth writing here. Its main 
physical features are to reveal that each term in the sum 
over n, accounts for the 4 possible transitions between the 
dressed states in the Jaynes-Cummings rungs n + 1 and 
n, as in the spontaneous emission case |12) . The first 
rung, or linear regime, is given by n = and consists 
of only the two transitions of the Rabi doublet. Other 
rungs give rise to a generalization of the Rabi frequency in 
the nonlinear regime, the nth-rung inner and outer Rabi 

frequencies: Ro,i[n] = y g 2 (V n +1± \fn) 2 — (P CT /4) 2 . 
In the Mollow triplet regime [P a > g), all the peaks 
positioned at the inner frequencies collapse at the cen- 
tre (including the Rabi doublet) giving rise to a single 
central peak. Outer peaks remain split at frequencies 
±R Q [n] « ±y/±g*n-\p a /4)*. 

The spectrum obtained with the previous derivations 
can be further simplified for the range of parameters in 
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FIG. 2: (Color online) (a) Comparison between our analyt- 
ical expression Eq. nty, without the elastic peak, (solid red) 
and the exact numerical solution (dashed blue), (b) Differ- 
ence between the incoherent (thick solid, red line) and the 
coherent (thin solid, black line) Mollow triplets, in equivalent 
conditions. The satellite peaks, that cause the departure, are 
plotted in dashed. Parameters: P„ — 6.3g and "/ a — O.lg. 



Eq. Q, to give a compact analytical expression. First, 
one considers only the coefficients with leading terms in 
n, making use again of n + 1 ~ n. Furthermore, due 
to the Poissonian statistics, only rungs with n close to 
n a contribute significantly to the spectra allowing the 
substitution n — > n a in Q[n\. The sum over n simplifies 
thanks to the normalization of the distribution function: 
Sn^N = 1- Finally, we neglect terms related to j a 
before those related to much larger rates, P a and K a , 
i.e., we write the spectrum for these three rates only, 
through the substitution g 2 = K CT 7 a /4, and then simply 
set 7 a — > 0. This results in the expression for Si nc (uj) in 
terms of P a and K a only: 
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This is our main result. The structure of the lineshape is 
the same as that of its coherent counterpart, Eq. (HI: a 
Dirac 8 function from the elastically scattered laser light 
superimposed to a triplet. The central Lorentzian peak 
has the same weight 1/2 but with FWHM given by the 
pump, P a . The two satellite peaks sit at ±#t(Ro) with 
Mollow splitting: 



Pn 



R Q = 1 fV8Ha/Pa-9 



(8) 



and FWHM 2>P a /2. The excellent agreement of our for- 
mula with exact numerical results jT2l [17] is shown in 
Fig. u[a) , where we superimpose in dashed blue the nu- 
merical computation to, in solid red, the analytical ex- 
pression Eq. ^h. Note the elastic peak in the numeric as 
a very narrow central line. 

Despite their similar structure, the lineshapes are in- 
trinsically of a different nature, as can be seen mathe- 
matically when reduced to their simplest, dimensionless 
expression, where they depend only on one parameter 
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FIG. 3: (Color online) Evolution of the Mollow triplet when 
increasing (from bottom to top) (a) the coherent excitation 
{Q.l/'Yct from 1 to 5.5 by steps of 0.5) and (b) the incoherent 
excitation {Pa lg from 2 to 11 by steps of 1, with -y a — O.lg). 



(£IlH(t and Po/ko)- Their similarities and differences 
can be better appreciated when compared on physical 
grounds, when the laser intensity for the conventional 
Mollow triplet, Eq. (fTl), is taken the same as the aver- 
age population of the cavity under incoherent excitation 
(i.e., we take fi| — > g 2 n a w Pa(n a — P a )/8), and when 
the broadening of the dot includes P a . In this way, we 
attempt the description of the incoherent system with 
the theory of the coherent one. We obtain an expression 
that shows the fundamental discrepancies between the 
two types of triplets: 
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Comparing this expression with Eq. d7J, one can see, 
indeed, that the central peak is the same in both cases, 
as well as the position and broadening of the satellite 
peaks (third terms have the same denominator), so the 
underlying structures bear much similarities. However 
the satellite lineshapes differ considerably, being strongly 
affected by the effect of incoherent pumping, that renor- 
malizes the rungs of the Jaynes-Cummings ladder, and 
other factors such as the dot population, which under 
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are shown in more details in Fig. ^[b), where they are 
plotted (in dashed) together with the whole triplets, in 
the coherent (thin black) and the incoherent (thick red) 
cases. 

Finally, Fig. [3] shows the natural experimental config- 
uration to demonstrate the new character of nonlinear 
spectroscopy in microcavities under incoherent pumping, 
and to contrast it with its coherent counterpart. Increas- 
ing pumping, one sees that in the coherent case (up- 
per panel), the triplet is better resolved, with a larger 
splitting, while in the incoherent case (lower panel), the 
satellites overlap with the central line as a result from 
pumping that splits them sublinearly, Eq. Q, and also 
increases their broadening. The two phenomenologies, 
despite their deep interconnections and common features, 
are strikingly different and the evidence of the new one 
should pose no problem even on qualitative grounds. 



[1 
[2 

[3 

[4: 

[5 
[6 

[7] 
[8] 

[9 

[io: 

in 

[12 
[13 

[14 

[15 
[16 



B. R. Mollow, Phys. Rev. 188, 1969 (1969). 
R. Loudon, The quantum theory of light (Oxford Science 
Publications, 2000), 3rd ed. 

F. Y. Wu, R. E. Grove, and S. Ezekiel, Phys. Rev. Lett. 
35, 1426 (1975). 

A. Muller, E. B. Flagg, P. Bianucci, X. Y. Wang, D. G. 
Deppe, W. Ma, J. Zhang, G. J. Salamo, M. Xiao, and 

C. K. Shih, Phys. Rev. Lett. 99, 187402 (2007). 
S. Ates, S. M. Ulrich, S. Reitzenstein, A. Lofflcr, 
A. Forchel, and P. Michler, Phys. Rev. Lett. 103, 167402 
(2009). 

E. B. Flagg, A. Muller, J. W. Robertson, S. Founta, D. G. 
Deppe, M. Xiao, W. Ma, G. J. Salamo, and C. K. Shih, 
Nat. Phys. 5, 203 (2009). 

A. N. Vamivakas, Y. Zhao, C.-Y. Lu, and M. AtatA^re, 
Nat. Phys. 5, 198 (2009). 

O. Astafiev, A. M. Zagoskin, J. A. A. Abdumalikov, Y. A. 
Pashkin, T. Yamamoto, K. Inomata, Y. Nakamura, and 
J. S. Tsai, Science 327, 840 (2010). 

S. Haroche and J.-M. Raimond, Exploring the Quantum: 
Atoms, Cavities, and Photons (Oxford University Press, 
2006). 

E. del Valle, Microcavity Quantum Electrodynamics 
(VDM Verlag, Saarbriicken, Germany, 2010). 
M. Nomura, N. Kumagai, S. Iwamoto, Y. Ota, and 
Y. Arakawa, Nat. Phys. 6, 279 (2010). 
E. del Valle, F. P. Laussy, and C. Tejedor, Phys. Rev. B 
79, 235326 (2009). 

S. Miinch, S. Reitzenstein, P. Franeck, A. Lofflcr, 
T. Heindel, S. Honing, L. Worschech, and A. Forchel, 
Opt. Express 17, 12821 (2009). 

A. Poddubny, M. Glazov, and N. Averkiev, 
arXiv: 1007.0244 (2010). 

Y. Mu and C. M. Savage, Phys. Rev. A 46, 5944 (1992). 
M. O. Scully and M. S. Zubairy, Quantum optics (Cam- 
bridge University Press, 2002). 

[17] M. Loffler, G. M. Meyer, and H. Walther, Phys. Rev. A 
55, 3923 (1997). 



